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A Numerical Study of Dedekind's Cubic Class Number 

Formula 

Harvey Cohn 1 

In this paper an analytic formula of Dedekind is used to compute class numbers for a 
sample of pure cubic fields, employing the National Bureau of Standards electronic com- 
puter, the SEAC. The computation is one of major magnitude, and it illustrates the useful- 
ness of combining integral and decimal arithmetic. The class numbers obtained are of 
much greater magnitude than those of pure cubic fields previously studied by means of 
hand computer methods. 



1. Introduction 

The discovery of fields with high class numbers 2 
has become an end in itself [10]. 3 There are gener- 
ally two possible procedures: One is the algebraic 
method, consisting, operationally, of a certain sorting 
of values of polynomials according to factors; the 
second is the analytic method, based on zeta-function 
residues. The second method is of special interest 
when applicable, because of the scarcity of types of 
field for which a usable formula emerges. 

In 1900, Dedekind applied both methods to the 
"pure" cubic field, or the field formed by adjoining 
M 1 ^ to the rationals [4]. The largest class number h 
that Dedekind achieved was h=9 for M=9l. In 
1950, Cassels published a table of class numbers [1] 
with the larger h=l2 for M=43, using the algebraic 
method. To achieve a substantially larger h it is 
necessary to take larger values of M. Furthermore, 
Dedekind's analytic formula is the simpler method 
to program for a machine for certain special M, 
namely, those that equal (s 3 +l)/c 3 for s and c 
integers. A sample of 21 cases were run, yielding 
six class numbers h between 18 and 27, as well as 
other smaller ones (including the cases cited above) . 

As an incidental result the quadratic fields (whose 
relation to the pure cubic field is similar to the rela- 
tion of rational numbers to the quadratic field) are 
studied in some numerical detail. 

The computation uses a mixture of integral and 
decimal arithmetic like the earlier work with Gorn 
[3], but, here, more fortunately, round-off provides 
no difficulty, even far beyond the range covered. 
The work was begun about July 1956 and the runs 
were completed November 1956. 

2. Dedekind's Formula 

The first algebraic investigations of class number 
were made by Markoff [8], but the first explicit 
"closed" formula for class number was derived by 
Dedekind. The entire numerical study undertaken 
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here stems from his formula, which we give without 
proof (in an equivalent form) : 



> = nH( % )^ 

6 HH(« S )W 



(2.1) 



Here h is the much desired class number of the 
real field R(MX), and e(>l) is the (often elusive) 
fundamental unit, which is discussed later. The 
function H(z) is defined for z a positive imaginary 
quantity (Im2>0), as the (positive) function 

H(2)=iKsM-s)tf(s-s)P*, (2.2) 

where rj(z) is the (analytic) Dedekind eta-function [6]. 



77(2) =exp {jiz\Vi) II (1 — exp 2Tizri), 

n = l 



(2.3) 



and z denotes the complex conjugate of z. The 
symbols u N , u R , are rather complicated and are 
treated in the next section, but more immediately 
relevant to the cubic field itself is the following fact: 
If M is decomposed by positive a, 6, c, so that 



M=ab 2 c 3 , 



(2.4) 



(so that a and b are free of quadratic factors), then 
R(M^) is generated equally well by (ab 2 )^ or 
(a 2 b)^. Hence the number pair (a,b) of relatively 
prime integers (not both 1) uniquely determines a 
field K a ,i,=R(MH). For this field the conductor k 
is defined as 



k-- 
k= 



-ab 
Sab 



As a final definition, 



D- 



-b 2 = (mod 9) 
-b 2 ^0 (mod 9) 

= -3£ 2 



.} 



(2.5) 



(2.6) 



is called the discriminant. Its significance lies in 
the fact that the algebraic integers, in particular the 
fundamental unit e, together with the algebraic 
conjugates e', 1' (the complex conjugates) satisfy 



£><|(«-e')(e'-I')(e-7')| S 



(2.7) 
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3. Quadratic Field 

Associated with the cubic field in question is 
the quadratic field generated by adjoining to the 
rationals the imaginary cube root of unity p given 
by (say) the positive imaginary root of 

p 2 +p+l = 0. (3.1) 

For this the norm function N is defined by 

N(p+q P ) = (p+<zp) (p+2P 2 ) =p 2 -pq+g 2 . (3.2) 

We finally associate with the field K atb and its 
conductor k the triples of integers {^4, B, C}, such 
that 

AC=k (3.3) 



1^B<A 

(A,N(B+C P ))-- 



= 1. 



(3.4) 
(3.5) 



For every such triple of integers we construct the 
quadratic form 



Q(x,y)=N[Ax+(B+C P )y] 



(3.6) 



and consider any prime p, not dividing 3k 2 , repre- 
sented by Q(x,y) for integral x and y. Then {A,B,C} 
is called a residue (or nonresidue) triple or Q(x,y) is 
called a residue (or nonresidue) form according to 
whether or not the equation 



x 3 =a 2 b (mod p) 



(3.7) 



is solvable for an integer x. Then for each triple 
satisfying eq (3.3), (3.4), or (3.5) a form Q(x,y) is 
defined with "root" 



B+Cp 



(3.8) 



The root co is an u R or an w N according as the residue 
or nonresidue character of the form. The product 
symbol in eq (2.1) should then be self-explanatory. 
It can be shown there are %k" triples {A,B,C} that 
generate an <a B and 6k" that generate an co N (not all 
of the o) R or co^ being distinct), where 9&", the total 
number of triples, is given by 



9k"=kIl(l-(-3/p)/p). 

Plk 



(3.9) 



Here (—3/p) is the usual quadratic residue symbol, 
1 for p = 1 (mod 3) 



(-3/*) = 



- 1 for p = — 1 (mod 3) 
0for^=0(mod3). 



(3.10) 



The properties of the quadratic form Q(x,y) will 
eventually come into play. It might suffice for the 



time being to mention that the function H(w), ac- 
cording to a well-known property of the Dedekind 
eta-function, is invariant under a unimodular change 
of variables; thus a unimodular integral change of 
variables of Q(x,y) will not affect H(co) ; neither will 
the change «'=•— c3. Thus let Q(x,y) be equivalent 
to the reduced form 



Q*(X,Y)=RX 2 +SXY+TY 2 . 
Its root may be written 

* -s+ihsy* 



(3.11) 



so that 



R 



0^S/2<R£(-Dl3)*=k 



(3.12) 



because, for all 9k" forms corresponding to con- 
ductor k, 



D=-4:RT+S 2 =-3k 2 . 



(3.13) 



We are now guaranteed by eq (3.12) a much more 
favorable scaling for the computation of H(w). 
Thus 



log H(«)= -Tk3^/(12R) + JJog(l-2K n cosSn/R+K 2n )\ 

K = exp(-TkS^/R) i n ^ ] 

(3.14) 

and k ^ exp( — ir3 1/2 ) , so that to the degree of round-off 
used here, it will develop that no more than five 
terms of the infinite sum (or product) are necessary, 
when the reduced form is used for Q(x,y). 

4. Units and Choice of Field 

We henceforth confine ourselves to fields generated 
by M 1/z , where, for a positive or negative integer s T 

M=s*+L (4.1) 

The advantage is that a very obvious unit is 

i 7 =|~8+Af 1/8 |- 1 = |« 2 +«Af 1/8 +M 2/8 | > (4.2) 

where M 1/3 refers to the real root. A fundamental 
unit is clearly of the type V /m > an( l the value of m 
can be ascertained in the sample treated here, by 
means of well-known results of Dedekind supple- 
mented by the lemma (compare [9]): 

Lemma. In the notation of (2.5), if &>10s, then 



e=r] or 7j 



j ii 



(4.3) 



Proof. Let e=r) 1/m be the fundamental unit. Then 
e, as an algebraic number, has a (root) discriminant 
not less than the discriminant D of K ttib . Hence 

D<\e-e'\ 2 
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where e', !' are the complex conjugate roots belong- 
ing to the same equation as e. But |e'| = |e'| = €- 1/2 , 
thus 

DS (6+€- 1/2 ) 4 (2e- 1/2 ) 2 ^4€ 3 (l + e- 3/2 ) 4 . (4.4) 

But finally, D=3k 2 , e=r] 1/m , whereas 77 is approx- 
imately 3s 2 from eq (4.2). Hence 

3F<4(3s 2 ) 3/M (l + r?- 3/2m ) 4 (V3s 2 ) 3/w . (4.5) 

Taking the maximum on the right when m^3, we 
see that the corresponding factors in parentheses 
are as follows: 



3F^4(3s 2 ).2 4 (l + 2 1/3 +2 2/3 )/3, 



(4.6) 



(the last factor being the one where m=3, s= 1), the 
inequality (4.6) is contradicted when fc>10s. Q.E.D. 
The list of 21 cases is given in table 1 in order of 
increasing k. Those cases in which the fundamental 
unit was known to Dedekind or Cassels are denoted 4 
by a "D" or "C", respectively. In the other cases, 
easily, k^>10s. The possibility that 7j=e 2 is excluded 
by the fact that in these cases the class number, 
calculated by assuming rj=e, will turn out to be odd. 

5. Subroutines 

The following number-theoretic subroutines were 
programed for this computation: One should bear 
in mind that the SEAC is a binary machine with a 
44 bit word and sign bit. The numbers treated 
range arithmetically from — 4+2~ 42 to 4 — 2~ 42 , i. e., 
the decimal point is fixed between the second and 
third high-order bit. Integers are held in the low 
end of the register (or scaled by 2 -42 ). 

1. Division of two integers to produce integral 
quotient and integral remainder. (The built-in 
division operation produced only a "decimal" 
quotient.) 

2. The Euclidean algorithm, or solution in X, Y of 



sX+tY=(s,t). 



(5.1) 



3. Reduction of modules. 5 

4. Multiplication of modules. 5 

5. Generation of coefficients of Q(x,y) (from { A,B,C} 
of eq (3.6)). 

6. Transformation of coefficients of Q(x,y) under 
unimodular change of variables. 

7. Reduction of Q(x,y) to the properly or impropeily 
eauivalent reduced form Q*(x,y) subject to conditions 
(3.12). To do this we first find the minimum value 
of Q(x,y) by scanning all of the finite set of points 
for which Q(x,yXQ(l,0)=A 2 , or for which 



[Ax+ (B+C/2)y] 2 +3C 2 y 2 /4;<:A 2 . 



(5.2) 



If the minimum value is Q(x ,y )=R, then the 
Euclidean algorithm generates an integral uni- 
modular transformation of coordinates, (or finds 



4 See headnote to table 1. 
* See section 7. 



#12,022), so that a 2 2X — a 12 y =l. Thus in the new 
coordinates, 



x X'+a l2 Y' = 

y,X'+a 22 Y' = 



f =y) 



(5.3) 



Q(x,y) becomes 



Q(X,Y)=RX /2 +S'X'Y' + T'Y' 2 . (5.4) 

Now to minimize the coefficient of X'Y 1 ', all we need 
do is perform the transformation 



X'=X+gY 1 



(5.5) 



choosing g=—[S'/2R], (the integral part), and we 
obtain Q*(X,Y)oi eq (3.11) and (3.12). 

8. Determination of the residue character of a Quad- 
ratic form. This subroutine finds a small 6 prime p 
represented by Q*(X,Y) and not dividing D. (The 
procedure is similar to the search for the minimum 
of Q(x,y) by eq (5.2).) Then the eq (3.7) is tested 
with -\-x and — x, where x runs from 1 to (p — 1)/2. 
If x passes through all these values without eq (3.7) 
being satisfied, then, and only then, Q is a non- 
residue form. 

9. The built-in u base zero" operations. These are 
standard fixed-point subroutines for input con- 
version to binary form and output conversion to 
decimal form, for logarithms, exponentials, sines, 
and cosines, among others. 

6. First Run. Class Numbers 

The subroutines listed earlier, taken together 
with the formulas and concomitant explanations, 
should indicate the basic program for computing h. 
It must be added that (binary) nonintegral quanti- 
ties used, e. g., in logarithmic computations were 
scaled by 2 -14 , and the final decimal values for 
h log 10 e were expressed to six decimal places. 

The basic procedure is to generate the 9fc" triples 
{A,B,C} in lexicographic order. For each triple we 
compute the corresponding forms Q and (reduced) 
Q*, and finally we form (and accumulate) the con- 
tributions (3.14), resulting in the following data (in 
decimal form) : 

Input: "Case number", a, b, k y D 

Output: For a nonresidue form, eq (3.11), 

R,S 



6 No effort was made to secure the smallest p, although this seems to have 
happened accidentally in the few cases spot checked. The smallest prime would 
have required much more programing, although the effect might have proved 
rewarding. 
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For a residue form, eq (3.7), (3.11), 

Final output, 

h log e, h logio e. 

For simplicity, R and S, which are never negative, 
were packed into one word, likewise x (with sign 
removed) and p. To identify a residue form, the 
x, p combination was printed first, with a minus 
sign. The residue and nonresidue forms are omitted 
from the final tabulation to save space. The 
division h=h logi €/]ogi e was performed by hand, 
resulting in integral h to five significant figures (see 
table 1). The data from tables of Dedekind and 
Cassels are denoted by "D" and "C", respectively. 

During the actual runs, the program was read in 
on wire, and the parameters were read in on teletype 
tape in trios of "case number", a, b. (Here the case 
number was packed with identifying information to 
check the tape position.) The output was on 
teletype printer. 

The program exhausted almost 900 of the 1,024 
electrostatic memory spaces. 

The running time varied from 3 minutes for case I 
to 14 hours for case XXI. A subroutine permitted 
the current {A,B,C} and cumulative sum of ±log 
H(co) to be dumped every half hour or so for rerun 
purposes. Almost no time went into input-output, 
by comparison. The vast bulk of the time went 
toward the testing of the residue characters, i. e., 
toward the solution or nonsolution of eq (3.7). The 
problem was brought to a close at case XXI largely 
because time was beginning to become prohibitive. 



The SEAC has no overflow stops at all, not even 
for 0/0, and no automatic error checks. As a substi- 
tute for both, a test was made to see if the discrimi- 
nant of each quadratic form checks with — 3k 2 . 
In the range covered, no overflow would have 
occurred, barring machine errors. 

7. Group Structure of Modules 

A good deal of insight into the machinery of 
Dedekind's formula is gained by considering the 
composition structure of the group of quadratic 
forms in the "order" of discriminant —3k 2 . Once 
again, rather than summarize a well-known and 
lengthy theory, we shall consider the modules only 
in relation to the present computation. 

The notation 

Wt=[A 1 +B lP ,...,A t +B t p] (7.1) 

denotes the aggregate of integers in R(p) 

f=Sn«(A,+B,p), 

where n t are integers. Such a set is called a module. 
Every module has a reduced basis like [A, B+Cp] 
identified with the triple {A,B,C} of section 3. 
Here A is taken as the minimum positive rational 
integer in 9JJ, and B+Cp is defined by the conditions 
that C is the minimum positive coefficient of p, 
whereas B is the minimum positive integer in a module 
element in which the coefficient of p is C. The enumer- 
ation of triples {A,B,C} enumerates modules, of 
course, whereas Q(x,y) is the norm of an arbitrary 
element of the module. Clearlv 



Table 1 
"D" and "C" denote units or class numbers appearing in Dedekind's and Cassel's works, respectively. 



Case 


s 


\\M\\ 


a 


b 


c 


logio v 


£ 




k 


k" 


h log 10 e 


h 




Gen- 
erator 


D 


I 


1 


2 


2 


1 


1 


0. 585159 


f\ 


D 


6 


1 


0. 585159 


1 


D 


21/ 3 


-108 


II 


2 


9 


1 


3 


1 


1. 096455 


V 


D 


9 


1 


1. 096455 


1 


D 


31/3 


-243 


III 


3 


28 


7 


2 


1 


1. 436650 


rji/2 


D 


14 


2 


2. 154975 


3 


D 


28 1 / 3 


-588 


IV 


-18 


5831 


17 


1 


7 


2. 987641 


V 


D 


17 


2 


2. 987641 


1 


D 


I71/3 


-867 


V 


8 


513 


19 


1 


3 


2. 283584 


V V2 


D 


19 


2 


3. 425376 


3 


D 


I91/3 


-1083 


VI 


-2 


7 


7 


1 


1 


1. 060137 


1? 


D 


21 


2 


3. 180412 


3 


D 


71/3 


-1323 


VII 


-3 


26 


26 


1 


1 


1. 425922 


V 


D 


26 


4 


4. 277770 


3 


C 


26 1 / 3 


-2028 


VIII 


19 


6860 


5 


2 


7 


3. 034650 


rji/2 


D 


30 


6 


4. 551974 


3 


D 


20 1 / 3 


-2700 


IX 


-10 


999 


37 


1 


:'» 


2. 476964 


f] 


D 


37 


4 


7. 430929 


3 


C 


371/3 


-4107 


X 


-4 


63 


7 


3 


1 


1. 678966 


T) 


D 


63 


6 


10. 073793 


6 


D 


63 1 / 3 


-11907 


XI 


-9 


728 


91 


1 


2 


2. 385408 


V 




91 


8 


21. 468668 


9 


D 


911/3 


-24843 


XII 


23 


12168 


1 


39 


2 


3. 200589 


V 


D 


117 


12 


19. 203533 


6 


C 


39 1 / 3 


-41067 


XIII 


5 


126 


14 


3 


1 


1. 876216 


V 




126 


18 


16. 885942 


9 




126i/ 3 


-47628 


XIV 


7 


344 


43 


1 


2 


2. 167739 


V 


D 


129 


14 


26. 012867 


12 


C 


431/3 


-49923 


XV 


-5 


124 


31 


2 


1 


1. 873900 


1) 




186 


30 


16. 865096 


9 




1241/3 


-103788 


XVI 


4 


65 


65 


1 


I 


1. 683490 


V 


D 


195 


24 


30. 302812 


18 




65i/3 


-114075 


XVII 


-6 


215 


215 


1 


i 


2. 032753 


V 




215 


28 


42. 687799 


21 




2151/3 


-138675 


XVIII 


6 


217 


217 


1 


1 


2. 034093 


t\ 




217 


20 


54. 920504 


27 




2171/3 


-141267 


XIX 


-7 


342 


38 


3 


I 


2. 166895 


f\ 




342 


54 


58. 506160 


27 




3421/3 


-350892 


XX 


15 


3376 


422 


1 


2 


2. 829347 


V 




422 


70 


59. 416280 


21 




4221/3 


-534252 


XXI 


17 


4914 


182 


1 


3 


2. 938049 


V 




546 


72 


79. 327311 


27 




1821/3 


-894348 
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P m=p[A,B+C P ] = [pA, -C+(B-C)p] 

is another module leading to the same Q(x,y), as 
N(p) = l, so is p'Wl. Thus the triples {A,B, C) or 
the reduced modules [A, i?+(7p], with AC=k, 
come in "threes," leading to the same form Q(x,y). 
Dedekind considered the multiplicative group 
formed by the modules under multiplication of basis 
elements: 

[A ) B 1 +C lP ][A 2j B 2 +C 2P ]= 
[A,A 2j A 1 B 2 +A l C 2 p 9 A 2 B 1 +A 2 C lPy (B 1 B 2 -C l C 2 ) + 
(B 1 C 2 +B 2 C l -C 1 C 2 )p]. (7.2) 

It can be shown that the product module can be 
reduced to the form [A z , B 3 +C z p] where A^C 2 =k, 
(provided A 1 C i =A 2 C 2 =k). This process is closely 
analogous to the composition of forms, or, of course 
the multiplication of ideals. 

For the purpose of performing this multiplication 
of modules we need two subroutines mentioned 
earlier (numbers 3 and 4 of section 5). 

3. Reduction of modules, (from [A u JBi + Cip, 
B 2 +C 2P ] to the reduced form [A, B+Cp]). 

4. Multiplication oj modules . This is the repeated 
application of the earlier subroutine for the reduction 
of product (7.2) to [A } B+Cp}. 

The notation "residue" (or "nonresidue") module 
can now be introduced to correspond to "residue" 
(or "nonresidue") form. Dedekind has shown that 
the set of residue modules forms a subgroup of index 
three among the set of all modules. Thus the non- 
residue modules can be partitioned into co-sets so 
that all modules are of type 

where 9? is a variable residue and % is a fixed non- 
residue module. In the notation 11(99?) for 
U(B-\-Cp/A), we can write in accordance with (2.1) 



IIHflfflO 



(7.3) 



Now if we take precautions that once SR appears, p9? 
and p 2 9? shall not appear, we have only k n factors in 
numerator or denominator. 

Now (7.3) is very close to Dedekind J s original 
formula. (Dedekind used quadratic forms and not 
modules to generate the products in (7.3).) 

8. Ambiguous Forms and Modules 

To appreciate the relation between modules and 
forms we must realize that not only do 9ft, p®^ p 2 9ft 
correspond to the same form but 9ft and 9ft can 
correspond to the same form. Here 9ft, the conju- 
gate module is defined by 



Wt=[A,B+Cp] 
m=[A,B+C P 2 ] 



1 



(8.1) 



The difficulty really arises when 9K=2K, (or P 2K=2tt, 
etc.). (The second possibility is really no different 
as it implies p 2 9ft is then its own conjugate.) These 
are the ambiguous forms of classical number theory. 
Generally 9ft9ft=[l, kp] the unit element, or principal 
ideal so that if 2W=W, 9ft 2 is a unit element and 9ft 
must belong to any subgroup of index three. In 
other words the self-conjugate modules are all 
residue modules. 

To give an ample illustration let us take case VI, 
s= — 4, k=21, and list the 9k" = 18 modules in lexi 
cographic order, as does the computing machine, but 
with the interrelations in question clearly indicated. 

List of modules for k=21 



[l f 21 P ] = 2Hb 
[3, 7 P ] = 3tt, 
[3, l+7p]-Mi 
[7, 3p]=pMi 



[7, 5 + 3p] = % 
[7, () + 3p] = W 2 
[21, p] = p9ft 



[21, 7 + p]=p9ftt 

[21,9 + p] = p9ft 3 
[21, 13 + P ]=p 2 ^ 



[21, l+p] = p 2 9ft [21, 15 + p]=p 2 9ft! 



[7, 3 + 3p] = p 2 2»i [21,4 + p] = p 2 9ft 2 [21, 16 + p]=p 2 SW 2 
[7, 4 + 3p] = Tl 2 [21, Q + p] = p¥l 2 [21, 18+p]= P 9ft 2 . 

These modules form a group with identity element 
9ft . Rather than construct the whole group table 
we might identify 9ftz,p9ft*,p 2 9ft* as one element, and 
it is easy to see that the group reduces to an abelian 
group of order 6 given by 

9%, 2Ri, %, Mi=2»?, m 2 =mm, W 2 =mm 2 (8.2) 

with the generators 9ft b 9ft 3 satisfying 

We see 9ft , 9ft 3 are the residue subgroup, for clearly, 

9fto=Ifto, 9ft3=9ft3; 

The corresponding quadratic forms, (grouped after 
Dedekind with their "negatives"), arc 7 

9ft : i?±xy+331y 2 =Q (x,y) 

9ft! : 9x 2 ±3xy+Z7y 2 =Q l (x,y) 

9ft 2 : 13x 2 ±9xy+27y 2 =Q 2 (x,y) 

9ft 3 : 19x 2 ±llxy+19y 2 =Q s (x,y).J 

Forms for 9ft and 9ft 3 are ambiguous (or self-conju- 
gate). Note that ab 2 =7, so that for the primes given 

f 331 = Qo(0,l), 10 3 = 7(mod 331) 

19 = ^(1,0), 4 3 =7(mod 19) 



(8.4) 



7 The negative of Q(x,y) is found by replacing S by —Sor, effectively, by keeping 
R and D fixed and replacing S by 2R-S. 
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Table 2 











Form 




Case 


Test value 










Residue 


Nonresidue 


\x\ 


7> 


R 


8 


IV 


S' 


I 


4 


31 


1 





4 


2 


II 


16 


61 


1 


1 


7 


3 


III 


! u . 


151 


1 





4 


2 


61 


3 





12 


18 


IV 


{ ,i 


223 


1 


. 1 


7 


13 


73 


3 


3 


13 


15 


V 


f 31 

I 26 


271 


1 


1 


7 


3 


97 


3 


3 


13 


23 


VI 


I '! 


331 


1 


1 


9 


15 


19 


a 19 


27 


13 


17 




f 49 


523 


1 





4 


6 


VII 


17 


181 


3 





12 


6 


4 


19 


19 


10 


21 


24 




I 4 


19 


19 


28 


7 


10 




( U 


691 


1 





4 


2 




22 


127 


25 





13 


2 


VIII 


- I 


19 
19 


19 
19 


6 
32 


25 

9 


20 
6 




1 


7 


7 


10 


25 


40 




I 1 


7 


7 


4 


27 


36 




f 304 


1, 033 


1 


1 


7 


11 


IX 


1 104 


349 


3 


3 


21 


39 


1 


19 


19 


23 


31 


27 




I 1 


19 


19 


15 


13 


25 




f 305 


3,019 


1 


1 


9 


15 




21 


73 


49 


49 


19 


33 


X 


< I 


37 


37 


9 


49 


63 


37 


37 


65 


43 


45 




1 


31 


31 


11 


49 


7 




I 1 


31 


31 


51 


13 


1 




( 747 


6, 211 


1 


1 


49 


7 




17 


139 


49 


49 


49 


77 




8 


67 


67 


125 


37 


13 


XI 


i 11 


79 


* 79 


147 


43 


21 


8 


67 


67 


9 


57 


3 




454 


2,089 


3 


3 


61 


17 




4 


31 


31 


53 


19 


3 




\ 4 


31 


31 


9 


73 


59 




3,364 


10, 267 


1 


1 


9 


15 




2 


139 


81 


63 


81 


9 




2 


139 


81 


99 


81 


45 




6 


193 


63 


39 


7 


11 




6 


193 


63 


87 


63 


3 


XII 


< I 5 


103 


» 103 


37 


49 


17 


1 


19 


19 


31 


67 


65 




3 


43 


43 


27 


79 


31 




4 


31 


31 


47 


73 


55 




4 


31 


31 


15 


37 


59 




3 


43 


43 


59 


61 


13 




1 


19 


19 


7 


97 


35 



a Signifies quadratic form whose negative belongs to a conjugate module. 









Form 






Test value 










Case 


Residue 


Nonresidue 


M 


V 


R 


8 


R' 


S' 




352 


11, 923 


1 





4 


2 




18 


331 


36 


66 


36 


30 




18 


331 


36 


6 


9 


6 




71 


439 


49 





73 


48 




43 


163 


81 


72 


79 


82 




43 


163 


81 


90 


43 


82 




2 


67 


67 


70 


81 


36 




2 


67 


67 


64 


49 


70 


XIII 


12 


103 


103 


24 


49 


84 


12 


103 


103 


182 


61 


108 




13 


109 


109 


130 


37 


18 




13 


109 


109 


88 


81 


18 




4 


31 


31 


22 


117 


102 




4 


31 


31 


40 


97 


22 




4 


19 


19 


10 


76 


10 




11 


241 


52 


54 


13 


2 




11 


241 


52 


50 


52 


102 




4 


19 


19 


28 


76 


66 




93 


25, 921 


1 





4 


6 




87 


751 


36 


66 


9 


12 




87 


751 


36 


6 


36 


42 




17 


1,069 


27 





108 


54 




109 


499 


52 


2 


13 


2 




55 


997 


63 


6 


103 


200 




64 


397 


76 


22 


76 


98 




56 


601 


148 


206 


148 


58 




56 


601 


148 


90 


37 


16 




64 


397 


76 


130 


19 


16 




55 


997 


63 


120 


163 


80 




109 


499 


52 


102 


52 


50 




3 


97 


97 


180 


157 


254 




12 


463 


117 


102 


61 


20 


XV 


3 


151 


151 


292 


43 


10 


7 


73 


73 


96 


91 


50 




163 


1,447 


91 


76 


79 


106 




13 


211 


133 


92 


109 


78 




31 


193 


171 


174 


49 


22 




23 


241 


108 


18 


27 


18 




115 


937 


28 


34 


7 


6 




115 


937 


28 


22 


28 


50 




23 


241 


108 


98 


108 


90 




31 


193 


171 


168 


133 


22 




13 


211 


133 


174 


63 


78 




163 


1,447 


91 


106 


117 


24 




7 


73 


73 


50 


139 


242 




3 


151 


151 


10 


172 


162 




12 


463 


117 


132 


127 


142 




3 


97 


97 


14 


67 


120 




' 838 


28, 549 


1 


1 


9 


3 




253 


1,297 


25 


25 


133 


191 




11 


211 


169 


169 


61 


103 




111 


409 


169 


273 


31 


37 




111 


409 


169 


65 


109 


131 




43 


157 


157 


165 


25 


35 




7 


139 


139 


105 


169 


143 




7 


139 


139 


173 


25 


5 




43 


157 


157 


149 


127 


239 




2 


19 


19 


1 


169 


299 




2 


19 


19 


37 


171 


. 75 


XVI 


6 


151 


151 


9 


169 


221 


18 


79 


79 


157 


49 


61 




18 


79 


79 


1 


169 


247 




6 


151 


151 


293 


175 


215 




12 


163 


163 


5 


37 


25 




2 


73 


73 


103 


175 


145 




48 


463 


63 


75 


63 


33 




4 


43 


43 


45 


103 


7 




342 


1,063 


27 


27 


27 


45 




4 


43 


43 


41 


133 


37 




48 


463 


63 


51 


7 


9 




2 


73 


73 


43 


67 


115 




I 12 


163 


163 


321 


97 


71 



270 



in accordance with eq (3.7). In the "general" situa- 
tion, of course, the residue modules will not all be 
ambiguous. 

9. Second Run. Residue Group 

A second run was undertaken to simulate more 
closely DedekincTs original method, which consisted 
of using the structure of the residue group. The 
advantage in so doing is clear if we think of the fact 
that most of the computing time goes toward testing 
eq (3.7) for a solution. The amount of work can be 
cut to about a sixth 8 if we agree to keep the quad- 
ratic forms in storage so as to never test the same 
form twice. Because the memory had been almost 
exhausted by the program, magnetic tape was used 
for storage of forms. 

The second run proceeded as follows. The ma- 
chine generated the modules, testing each in turn 
until it came to (say) %, the first nonresidue module, 
which it put into memory. (Only the internal 
memory had been used thus far.) Once the machine 
reached this point it returned to the beginning of the 
list of modules and did two things "simultaneously": 
First of all, it tested each module whose quadratic 
form was not on the tape and loaded on the tape 9 
every form that had been tested. Second, the ma- 
chine took the residue module (say) 9t that had not 
been previously tested (and loaded on tape) and 
formed fti?(t loading both quadratic forms on tape, 
while it accumulated log H(5fti9i)— log H(5R). 

The second run led to an interesting array of 
quadratic forms which we see in table 2. Here R,S 
refer to coefficients of the residue quadratic form 
and R',S' refer to the corresponding nonresidue 
quadratic form which is formed by multiplying by 
the module 5fti as explained earlier. 

Unfortunately, here overflow was a serious prob- 
lem in the multiplication of modules, and only cases 
I through XVI, excluding XIV, went to completion. 
The forms for which S=R or zero obviously correspond 
to self-conjugate modules. 10 The other self-conju- 
gate modules are marked by " a ", and can be spotted 
by the fact that the residue quadratic form and its 
negative do not both appear in table 2 (see footnote 
9). 

8 Unfortunately tape motion produced too much delay (as much as 1 min per 
search) to permit these benefits to be fully realized. 

9 Here, whenever 2ft = 9ft, p2ft, or p 2 9tt, one should load both the for 9tt and the 
negative (for M) on the tape to prevent both Q and its negative from appearing 
as residue forms. This fact had been overlooked and had to be corrected by an 
undesirable manual intervention in the computation. 

i° See footnote to table 2. 



10. Concluding Prospects 

The computation consummated here represents 
one phase of what really amounts to a larger task. 
The computation of units should also be programed, 
possibly after Voronoi's algorithm [5], so that h 
can be found even when M does not have the for- 
tuitous form s 3 +l. The storage of the previously 
tested forms in a larger high-speed memory would, 
of course, have speeded up computations, but the 
testing of congruence (3.7) leaves much to be desired 
in the direction of efficiency. To carry the groupings 
of modules and forms on to larger k, a double pre- 
cision multiplication of modules would have to be 
provided. Lastly, it would be desirable to find the 
class structure algebraically [2], which is not alto- 
gether beyond the capacity of the SEAC. 

The class structures, as well as the group struc- 
tures, of the modules shall be the subject of further 
theoretical investigations. 



The author thanks Olga Taussky for her aid and 
encouragement. 
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